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Abstract 

A diagrammatic method is presented for averaging over the circular ensemble 
of random-matrix theory. The method is applied to phase-coherent conduc- 
tion through a chaotic cavity (a "quantum dot") and through the interface 
between a normal metal and a superconductor. 
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I. INTRODUCTION 



The random-matrix theory of quantum transport describes the statistics of transport 
properties of phase- coherent (mesoscopic) systems in terms of the statistics of random ma- 
trices (for reviews, see Refs. 0-fl])- There exist two separate (but equivalent) approaches: 
Either the random matrix is used to model the Hamiltonian of the closed system, or it is 
used to model the scattering matrix of the open system. The second approach is more direct 
than the first, because the scattering matrix directly determines the conductance through 
the Landauer formula, 

G = ^trttl (1.1) 
h 

(The transmission matrix t is a sub matrix of the scattering matrix.) 

Random-matrix theory has been applied successfully to two types of mesoscopic systems: 
chaotic cavities and disordered wires. Baranger and Mello || and Jalabert, Pichard, and 
Beenakker || studied conduction through a chaotic cavity on the assumption that the scat- 
tering matrix S is uniformly distributed in the unitary group, restricted only by symmetry. 
This is the circular ensemble, introduced by Dyson |7[], and shown to apply to a chaotic 
cavity by Bliimel and Smilansky || . The symmetry restriction is that SS* = 1 in the pres- 
ence of time-reversal symmetry. (The superscript * indicates complex conjugation if the 
elements of S are complex numbers; in the presence of spin-orbit scattering, S is a matrix 
of quaternions, and S* denotes the quaternion complex conjugate.) For the disordered wire, 
the circular ensemble applies not to the scattering matrix itself, but to the unitary matrices 
v, w, v', and w' in the polar decomposition, 




The matrix T is a diagonal matrix containing the transmission eigenvalues T n e [0, 1] on 
the diagonal. (The T n 's are the eigenvalues of the matrix product tv .) The distribution 
of the transmission eigenvalues is governed by a Fokker-Planck equation, the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation PJIUf. The isotropy assumption [IIJ states that v, 



v', w, and w' are uniformly and independently distributed in the unitary group, with the 
restriction v*v' = 1, w*w' = 1 in the presence of time- reversal symmetry. 

The role of the circular ensemble of unitary matrices in the scattering matrix approach is 
comparable to the role of the Gaussian ensemble of Hermitian matrices in the Hamiltonian 
approach. However, whereas many computational techniques have been developed for aver- 
aging over the Gaussian ensemble [|TlH T8| , the circular ensemble has received less attention. 
If the dimension N of the unitary matrices is small, the average over the circular ensemble 



can be done exactly ||19| , p0| . For some applications in the regime of large N, one may regard 



the elements of the unitary matrix as independent Gaussian variables and then use 



the known diagrammatic perturbation theory for the Gaussian ensemble []12| , |T7|| . In other 
applications the Gaussian approximation breaks down. 

In this paper we present a diagrammatic technique for integration over the unitary group, 
which is not restricted to the Gaussian approximation. We discuss two applications: A 
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chaotic cavity coupled to the outside via a tunnel barrier, and a disordered wire attached to a 
superconductor. In both cases, we calculate the mean and variance of the conductance up to 
and including terms of order 1. We point out the analogy between the diagrams contributing 
to the average over the circular ensemble and the diffuson and cooperon diagrams which 
appear in the theory of weak localization 



22.231 and universal conductance fluctuations 



24] , |25|| in disordered metals. In the presence of the superconductor a third type of diagrams 



shows up, which gives rise to the coexistence of weak localization with a magnetic field 
2~6| , |27|1 , and to anomalous conductance fluctuations 



The paper starts in Sec. [TT] with a summary of known results [5U-fm for the integration 
over the unitary group of a polynomial function of matrix elements. The diagrammatic tech- 
nique is explained in Sec. JT|. Generalizations to unitary symmetric matrices and to unitary 
quaternion matrices are given in Sees. [TV] and [V], respectively. We then apply the technique 
to the chaotic cavity (Sec. |VT|) and the normal-metal-superconductor junction (Sec. [VI H ). 
Some of the results of Sec. |V| have been obtained previously by Iida, Weidenmuller, and 
Zuk, who used the Hamiltonian approach to quantum transport and the supersymmetry 
technique [1],|15|]- The resu lts of Sec. |V11| have been published in Refs. p6| , [28[] , without the 
detailed derivation presented here. There is some overlap between Sec. |VII| and a recent 
work by Argaman and Zee 



II. INTEGRATION OF POLYNOMIALS OF UNITARY MATRICES 



In this section we summarize known results p0|-p2| for the integration of a polynomial 
function f(U) of the matrix elements of an iV x N unitary matrix U over the unitary group 
U(N). We refer to the integration as an "average", which we denote by brackets (•■■), 

(/> = JdUf(U). (2.1) 

Here dll is the invariant measure (Haar measure) on U(N), normalized to unity (/ dU = 1). 
The ensemble of unitary matrices that corresponds to this average is known as the circular 
unitary ensemble (CUE) [0,[£J. 

We consider a polynomial function f(U) = U aih . . . U anbn U* lPl . . . U* mf3m . The average 
(/({/)) is zero unless n = m, ai, . . . , a n is a permutation P of a,\, . . . , a n , and (3i, . . . , /3 n is 
a permutation P' of b%, . . . , b n . The general structure of the average is 



U aibi ■ ■ ■ U a m bJ J l ll p 1 ■ ■ ■ U Lf3n) ~ 5nm 12 Vp . p ' IT ^'"PO-A^P'OV ( 2 - 2 ) 

P,P> j=l 

where the summation is over all permutations P and P' of the numbers 1, ...,n. The 
coefficients Vp,p/ depend only on the cycle structure of the permutation P~ 1 P' |^T|. Recall 
that each permutation of 1, . . . , n has a unique factorization in disjoint cyclic permutations 
("cycles") of lengths ci, . . . , Ck (where n = J2j=i c k)- The statement that Vp,p' depends only 
on the cycle structure of P~ 1 P' means that Vp^pi depends only on the lengths ci, . . . , c& of 
the cycles in the factorization of P~ 1 P'. One may therefore write V Cli ... jCfc instead of Vp t p/. 

As an example, we consider the case n = m = 2 explicitly. The summation over the 
permutations P and P' extends over the identity permutation id = [(1, 2) — >• (1, 2)] and the 
exchange permutation ex = [(1,2) — > (2, 1)]. Hence Eq. (|2.2|) reads 
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\^ai&i^a 2 &2^ai/3i ^22/32/ ~~ Vid,id ^aiai^biPi ^0202 ^f>2/32 + ^aia2^6i/3i^a 2 a!i^&2/32 

4" Vid,ex ^aiai^6i/32^a2a2^2/3i ^ex,ex <5(iia2 ^6i/32^«2Ql ^&2/9l ■ (^'^) 

The permutation P~ 1 P' that corresponds to P = P' = id [the first term on the r.h.s. of Eq. 



is again the identity permutation: P~ X P' = id = [(1,2) — > (1,2)]. Its factorization 
in cyclic permutations is id — (1 — > 1)(2 — > 2), so that P~ l P' factorizes in two cyclic 
permutations of unit length. Hence the cycle structure of P~ l P' is {1,1}, and V^id = 
14,1. The second term on the r.h.s. of Eq. ( |2.3|) , corresponding to P = ex, P' = id, has 
P~ 1 P' = ex = [(1, 2) — > (2, 1)], which factorizes in a single cyclic permutation of length two, 
ex = (1 — > 2 — > 1). Hence the cycle structure of P~~ X P' is {2}, and V^id = V 2 . Treating the 
remaining two terms of Eq. ( [2.3|) similarly, we obtain 



^1161^262^1^1^2/32/ _ ^1.1 ^aioi^ift^aaaa^afe + ^2 <W2<Wl^a 2 £*i<W2 

^aiai^6i/32^a 2 a2^62/3i + ^1,1 ^ai^^i/fe^ai^Wi • (2-4) 

In general, the coefficient Vi 1 refers to equal permutations P = P' , corresponding to a 
pairwise (Gaussian) contraction of the matrices U and U*. Coefficients V Clr .., Ck with some 
Cj 7^ 1 give non-Gaussian contributions. 



The coefficients V are determined by the recursion relation [[31 



k 

NV C \,...,Cff + ^ Vp,q,C2, — ,C k ~l~ C ,7 ^Cl+C,-.C 2 . — ,C,--1,Cj-|-1,...,Cfe ~ b^lVci Cfc) (2-5) 

P+9=ci j'=2 

with Vo = 1. One can show that the solution V cl ^ Ck does not depend on the order of the 
indices c\, . . . ,Ck- Results for V up to n = 5 are given in App. [A[ The large- iV expansion of 
V is 

Ki,..., Cft = n^. + C(iV fc - 2ri - 2 ), (2.6a) 
i=i 

V c = l -N l - 2c {-lf- 1 + 0{N- l ~ 2c ). (2.6b) 

(The numbers c -1 ^ 26 ^ 2 ) are the Catalan numbers.) For example, Vi y .. ) % = N~ n + (D(N~ n ~ 2 ). 
The Gaussian approximation amounts to setting all V^'s equal to zero except V\ y „^, which 
is set to N~ n . 

The coefficients V Cu ... )Ck determine the moments of U. Similarly, the coefficients W Cu ... )Ck 
determine the cumulants of U. The cumulants are obtained from the moments by subsequent 
subtraction of all possible factorizations in cumulants of lower degree. For example, 

W C1 = V C1 , (2.7a) 

^01,02 = ^1,02-^01^, (2.7b) 

^01,02,03 = Ki,o 2 ,o3 - W Cl W C2tC3 - W C2 W cliC3 - W C3 W CUC2 - W C1 W C2 W C3 . (2.7c) 
The recursion relation (|2.5[) for V implies a recursion relation for W, 
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k 

NW C i,...,c k + Wp,g,c 2 ,...,c fe + ^ Cj W C i+Cj,C2,—,Cj-i,Cj +1 ,...,ci t 

P+9=ci j=2 

k 1 

+ X] 7r TTTTT nt X! ^ / p.cp(2)v,c P(i) W / (?,cp (;+1) ,c P(fe) = 0, (2.8) 

p+q=ci Z=l 1/ _ J-J 1 ^ _ i J ! P 

with VFq = 1 and P a permutation of 2, . . . , k. To leading order in 1/N this equation has 
the solution, 

W Cu ..., Ck = 2^- 2 "- fc+2 (-l)" +fc (2n + fc ~ 3)! fl + 0(iV" 2 "- fc ). (2.9) 

Notice that W clt _ >Ck decreases with increasing number of cycles k, opposite to the behavior 

° f V ci,...,c k - 

In principle, the recursion relations permit an exact evaluation of the average of any 
polynomial function of U. In practice, as the number of IPs and U* : s increases, keeping 
track of the indices and of the Kronecker delta's which connect them becomes more and 
more cumbersome. It is by the introduction of a diagrammatic technique that one can carry 
out this bookkeeping problem in a controlled and systematic way. 



III. DIAGRAMMATIC TECHNIQUE 

The usefulness of diagrams for the bookkeeping problem is well-established for averages 
over the Gaussian ensemble of Hermitian matrices |12| . Brezin and Zee [17] have devel- 



oped a diagrammatic method which can be applied to non-Gaussian ensembles as well, as 
a perturbation expansion in a small parameter multiplying the non-Gaussian terms in the 
distribution. No such small parameter exists for the circular ensemble. The method pre- 
sented here deals with non-Gaussian contributions to all orders. Creutz |30| has given a 
diagrammatic algorithm for integrals over SU(iV). Because of the more complicated struc- 
ture of SU(iV), we could not effectively apply his method to integrals over U(N) in the case 
of a large number of £/'s. 

The diagrams consist of the building blocks shown in Fig. [1|. We represent matrix 
elements U ab or U** by thick dotted lines. The first index (a or a) is a black dot, the 
second index (6 or (3) a white dot. A fixed matrix Aij is represented by a directed thick solid 
line, pointed from the first to the second index. Summation over an index is indicated by 
attachment of the solid line to a dot. As an example, the functions f{U) = tr AUBW and 
g(U) = tr AUBUCUWW are represented in Fig. |. 

The average over the matrix U consists of summing over all permutations P and P' in 
Eq. ( |2.2| ) . Permutations are generated by drawing thin lines (representing Kronecker deltas) 
between all black dots attached to U and black dots attached to U* (one line per dot). 
Black dots connect to black dots and white dots to white dots. To find the contribution of 
the permutations P and P' to (f(U)), we need (i) to determine the cycle structure of the 
permutation P _1 P', and (ii) to sum over the indices of the fixed matrices A. 

(i) The cycle structure can be read off from the diagrams. A cycle of the permutation 
p-ip/ gj ves r j se to a closed circuit in the diagram consisting of alternating dotted and thin 
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FIG. 1. Substitution rules for the unitary matrices U and U*, the fixed matrix A and the 
Kronecker delta. 




FIG. 2. Diagrammatic representation of the functions f(U) = tr AUBU^ and 

g(U) = tvAUBUCU^DU^. 

lines. The length of the cycle is half the number of dotted lines contained in the circuit. 
We call such circuits [/-cycles of length c^. 

(ii) The trace over the elements of A is done by inspection of the closed circuits in 
the diagram which consist of alternating thick and thin lines. We call such circuits T- 
cycles. A T-cycle containing the matrices A^ , A^ , . . . , A^ (in this order) gives rise to 
trA^A^ . . . A( k \ If the thick line corresponding to a matrix A is traversed opposite to its 
direction, the matrix should be replaced by its transpose A T . 

As an example, let us consider the average of the functions f(U) = tvAUBW and 
g(U) = trAUBUCUWW. Connecting the dots by thin lines, we arrive at the diagrams of 
Fig. ^. For /, there is only one diagram. It contains a single [/-cycle of length 1 (weight 
Vi) and two T-cycles (which generate tr A and trB). We look up the value of V\ = 1/N in 
App. 0, and find 

(f(U)} = V 1 trAtrB = N-hrAtrB, (3.1) 

Four diagrams contribute to g. The first diagram contains two [/-cycles of length 1, and three 
T-cycles. Its contribution is Vi } itr Atr BD tr C . The second diagram contains two [/-cycles 
of length 1 and a single T-cycle. Its contribution is Vi^trADCB. The third and fourth 



6 




FIG. 3. Diagrammatic representation of the averages of the functions / and g in Fig. ^. 

diagram each contain a single [/-cycle of length 2 and two T-cycles. Their contributions are 
V 2 tr Ati BDC and V 2 tr ADB tr C. In total we find 

( g (U)) = V 1A (tiAtrBDtrC + tr ADCB ) + V 2 (tr A tr BDC +tvADBtvC) 
= (N 2 - l)' 1 (tr A tr BD tr C + tr ADCB) 

- [iV(iV 2 - l)]" 1 (trAtr J BDC , + trAD J BtrC). (3.2) 

Whereas each individual T-cycle gives rise to a trace of matrices, it is only the combi- 
nation of all [/-cycles together that determines the coefficient V C1 Ck . The evaluation of a 
diagram would be more efficient, if we could attribute a weight to an individual [/-cycle. We 
introduced the cumulant expansion of the coefficients V in the coefficients W for this pur- 
pose. The leading term V Cu ... jCk = ITp=i Wc p of the cumulant expansion attributes a weight 
W c to each individual [/-cycle of length c p . This is sufficient for the calculation of the 
large- N limit of the average (/). The next term Y<i<jW CitCj Ylp-ajW^ attributes a weight 
W CitCj to the pair of [/-cycles, and the weight W c to all others individually. This is 
sufficient for the variance of /. The general rule is that the j'th order cumulant of / in the 
large-iV limit requires the jth order term in the cumulant expansion of the coefficients V, 
and hence requires consideration of groups of j [/-cycles. 

Let us summarize the diagrammatic rules: 

1. Draw the diagrams according to the substitution rules of Fig. [I]. 

2. Draw thin lines to pair black dots attached to U to black dots attached to U*. Do the 
same for the white dots. 
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<f 2 )= 



v ■ ■ ■ \J B 1 " " Cr b 




B ■ ■ ■ (_7 B ' ' ' 

FIG. 4. Diagrammatic representation of (/■* 




3. Every closed circuit of alternating thick solid lines and thin solid lines (a T-cycle) 
corresponds to a trace of the matrices A appearing in the circuit. If a thick line is 
traversed opposite to its direction, the transpose of the matrix appears in the trace. 

4. Every closed circuit of alternating dotted and thin solid lines (a [/-cycle) corresponds 
to a cycle of length equal to half the number of dotted lines. The set of [/-cycles 
in a diagram defines the coefficient V C1 Ch , which is the weight of the diagram. The 
coefficient V can be factorized into cumulants. To determine the cumulant coefficients 
W, partition the [/-cycles into groups. Every group of p [/-cycles of lengths ci, . . . , c p 
contributes a weight W Clr „ jCp . 

The diagrammatic rules are exact. In the large- N limit, we may reduce the number of 
diagrams and partitions that is involved. Let us determine the order in N of a diagram 
with / T-cycles and k [/-cycles of total length n partitioned into g groups. Counting every 
trace as an order N and using the large-iV result ( j2.9| ) for the coefficients W, we find a 
contribution of order ]\r 2 g+/-A;-2n gj nce g < k the order is maximal if g = k and the total 
number of cycles k + 1 is maximal. Thus, for large N, we may restrict ourselves to diagrams 
with as many cycles as possible and with a partition of the [/-cycles in groups of a single 
cycle (i.e. we may approximate V clt ^ Ck ~ W C1 . . . W Ck ). 

We conclude this section with one more example, which is the calculation of the variance 
var/ = (f 2 ) — (f) 2 of the function f(U) = tr AUBU\ Diagrammatically, we calculate (f 2 ) 
as in Fig. [|a, resulting in 



+ W 2 



(f 2 ) = V 1A [{tr A) 2 (tr B) 2 + trA 2 tiB 2 
var/ = Wi,i [(tr A) 2 (tr B) 2 + tr A 2 tr B 2 + W 2 tr A 2 tr B 2 
+ W 2 [trA 2 (trB) 2 + (tr A) 2 tr B 



trA 2 (trB) 2 + (tr A) 2 tr B' 



(3.3a) 
(3.3b) 



If we now consider the order in iV of the various contributions, we see that the leading 0(N 2 ) 
term of (f 2 ) (I = 4, g = k = 2, corresponding to 6 cycles and a partition of the [/-cycles 
into two groups of a single cycle), is exactly canceled by (f) 2 - This exact cancelation is 
possible because the leading contribution of (f 2 ) is disconnected: Each T-cycle, and each 
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group of [/-cycles belongs entirely to one of the two factors trAUBW of f 2 . Only connected 
diagrams contribute to the variance of /. The connected diagrams are of order 1 (k + / = 4 
and g = k or k + I = 6 and g = k — 1). They give the variance 



var / = W 1A (tr A) 2 (tr B) 2 + W?tr A 2 tr B 2 

+ W 2 [tr A 2 (tr 5) 2 + (tr A) 2 tr 5 2 1 + ©(iV -1 ; 



(3.4) 



IV. INTEGRATION OF UNITARY SYMMETRIC MATRICES 

In the presence of time-reversal symmetry the scattering matrix S is both unitary and 
symmetric: SS^ = 1, S — S T . The elements of S are complex numbers. (The case of a 
quaternion S, corresponding to spin-orbit scattering, is treated in the next section.) The 
ensemble of uniformly distributed unitary symmetric matrices is known as the circular or- 
thogonal ensemble (COE) |7|,|33|. Averages of the unitary symmetric matrix U over the 
COE can be computed in two ways. One way is to substitute U = VV T , with the matrix 
V uniformly distributed over the unitary group. This has the advantage that one can use 
the same formulas as for averages over the CUE, but the disadvantage that the number of 
unitary matrices is doubled. A more efficient way is to use specific formulas for the COE, 
as we now discuss. 

The average of a polynomial in U and U* over the COE has the general structure 

2n 

P j=l 

The summation is over permutations P of the numbers 1, . . . , 2n. We can decompose P as 

p = I IT tA p c p ( ft rA , (4.2) 




where T 3 - and T- permute the numbers 2j — 1 and 2j, and P c (P Q ) permutes n even (odd) 
numbers. Because U a b = Ub a , the moment coefficient Vp depends only on the cycle structure 
{ci, . . . , Cfc} of P^Pq 01, so that we may write V^...^ instead of Vp. 



The moment coefficients obey the recursion relation 

k 

(N + Ci)V cl) ... jCjfe + ^ Vp,q,C2,...,c k + 2 Cj Ki+Cj,C2,...,Cj_i,Cj + i,...,c fc = ^cilK 2 ,...,c fc j (4-3) 
p+g=ci 3=2 

with Vo = 1. The large- expansion of V is 

V Cl _ Ck = f[V Cj +0(N k - 2n - 2 ), (4.4a) 
K = ijV 1 - 2 ^-!) - 1 ~ x 2 J - AA- 2c (-4) c - 1 + ^(Ar- 1 -^). (4.4b) 
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FIG. 5. Diagrammatic representation of (f(U)) for f(U) = ti AUBU\ where U is a unitary 
symmetric matrix. The second term arises because of the symmetry constraint. 



Compared with Eq. ( [2.6|) an extra term of order N 2c appears in V c because of the symmetry 
restriction. The recursion relation for the cumulant coefficients W is 

k 

(N + Cl )W Cu ..., Ch + £ W m ,c s ,...,c k + 2 Cj w / cl+CjiC2i ... iCj _ ljCj . +lv .. iCfc + 

p+q=c 1 j=2 
k x 

+ ^2 ^2jj~Z iy.(k — l)\ ^ ^ / P' C ^(2)v,Cp( i )^ / 9:Cp (i+1) ,Cp (fe) = 0, (4.5) 



with Wq = 1 and P a permutation of the numbers 2, . . . , k. The solution for large N is 

_ k 

r>2fc-l ftf-2n-k+2( _-( ^n+A; l^"- 1" ^ '-'J 1 

ci,...,c fe 



M = 2 2fc - 1 Ar-^ + 2 ( _ 1) n +fc (2n + fe 3)! ^ (2c, 1)1 + G(Ar _ 2ri _ fc+1) _ (4 g) 



3=1 



The coefficients V clr __ }Ck and W Clt ,„^ k are listed in App. [A] for n = c\ + . . . + Ck < 5. 

For the diagrammatic representation, we again use the substitution rules of Fig. |T[ The 
symmetry of U is taken into account by allowing thin lines between black and white dots. 
Therefore, rule 2 is replaced by 

2. Pair the dots attached to U to the dots attached to U* by connecting them with thin 
lines. 

As examples, we compute the averages of f(U) = trAUBW and g(U) = 
tr AUBUCW DW over the COE. The diagrams for (f(U)) are shown in Fig. |^, with the 
result 

(f(U)) = V^tvAtrB + trA T B) = (N + l) _1 (tr Ati B + tr A T B). (4.7a) 

Similarly, we find that 

(g(U)) = [(N + l)(N + 3)]- 1 

x (trAtrBDtrC + tr AD T B T tr C +trAtrBC T D + tr AD T CB T 
+ tr ADCB + tr AC T D T B + tr ADB T C T + tr AC T tr BD T ) 
- [(NiN + ^iN + S)]- 1 

x (tr A trB DC +trAC T D T B T +tiAtxBD T C + tr AC T DB T 

+ tr ADB tr C + tr AD T B tr C + tr ADBC T + tr AD T BC T 

+ ti AD T B T C T +ti AC T trBD + trAD T CB + trAC T DB 

+ trAtrBC' T D T +tr ADC B T + tvAtxBD T trC + tr ADB T trC). (4.7b) 
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V. INTEGRATION OF MATRICES OF QUATERNIONS 



We extend the results of the previous sections for integrals over unitary matrices of 
complex numbers to integrals over unitary matrices of quaternions. This is relevant to the 
case that spin-rotation symmetry is broken by spin-orbit scattering. 



Let us first recall the definition and basic properties of quaternions [ 33| . A quaternion q 
is represented by a 2 x 2 matrix, 

q — a l + %a\(j\ + ia 2 a 2 + ia 3 a 3 , (5.1) 

where II is the 2x2 unit matrix and cr } - is a Pauli matrix, 






The coefficients aj are complex numbers. The complex conjugate q* and Hermitian conjugate 
of a quaternion q are defined as 

q* = OqI + ia[a\ + %a^(J2 + ia^a^, q^ = a^TL — ia\ai — ia^o-i — ia^a^. (5.3) 

The complex conjugate of a quaternion differs from the complex conjugate of a 2 x 2 matrix, 
whereas the Hermitian conjugate equals the Hermitian conjugate of a 2 x 2 matrix. Let Q 
be an N x N matrix of quaternions with elements Qki = Qjy 1 + iQki a i + ^Qki a 2 + iQ^i a ^- 
The complex conjugate Q* and Hermitian conjugate are defined by (Q*)m = Qti an d 
(Q ] )ki = Q\ k - The dual matrix Q R is defined by Q R = (Qt)* = (Q*)t. We call Q unitary if 
QQ^ = 1 and self-dual if Q = Q R . A unitary self-dual matrix is defined by QQ^ = QQ* = 1. 
The trace trQ is defined by trQ = 52jQjj, which equals 1/2 the trace of the 2N x 2N 
complex matrix corresponding to Q. The scattering matrix in zero magnetic field is a 
unitary self-dual matrix, because of time-reversal symmetry. The ensemble of quaternion 
matrices which is uniformly distributed over the unitary group is called the circular unitary 
ensemble (CUE). If the ensemble is restricted to self-dual matrices it is called the circular 
symplectic ensemble (CSE) [0,[33|. 

The integration of a polynomial function f(U) of an N x N quaternion matrix U over 
the CUE or CSE can be related to the integration of a function f(U) of an N x N complex 
matrix U over the CUE or COE. The translation rule is as follows (a similar rule has been 
formulated for Gaussian ensembles in Refs. p5| , p6| ): 



1. f(U) is constructed from f(U) by replacing, respectively, the complex conjugates, Her- 
mitian conjugates, and duals of quaternion matrices by complex conjugates, Hermitian 
conjugates, and transposes of complex matrices. Furthermore, every trace is replaced 
by — |tr, and numerical factors iV are replaced by —\N. 

2. The average (f(U)) is calculated using the rules for integration of N x N complex 
matrices over the CUE or COE. 

3. The average (f(U)) over the CUE or CSE is found by replacing, respectively, the 
complex conjugates, Hermitian conjugates, and transposes of complex matrices by the 
complex conjugates, Hermitian conjugates, and duals of quaternion matrices. Traces 
are replaced by —2 tr and numerical factors iV by —2N. 
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As examples, we compute the averages of the functions f(U) = trAUBW and g(U) = 
tr AUBUCW DW of N x iV quaternion matrices over the CUE and CSE. The first step is 
to construct the functions f(U) and g(U) of N x N complex matrices, 

f (U) = -\ti AUBU\ g(U) = -\tiAUBUCU ] DU ] . (5.4) 
The second step is to average / and g over the CUE. The result is in Eqs. ( |3.1| ) and 



(/)cue = -\N~Ht AtxB, (5.5a) 
(g) CVE = -\(N 2 - I)' 1 (tiAtiBDtiC + tr ADCB) 

+ \[N(N 2 - l)]" 1 (tr A tr BDC + tr ADB tr C) . (5.5b) 

The third step is to translate back to quaternion matrices, 

(/)cue = N-hiAtiB, (5.6a) 
(g) CVE = (AN 2 - I)" 1 (4 tr A tr BD tr C + tr ADCB) 

- [N(AN 2 - I)]" 1 (tiAtiBDC + tiADBtiC) . (5.6b) 

Similarly, to average of / and g over the CSE we need the average of / and g over the COE 
given by Eq. ( f4.7a|) , and then translate back to quaternion matrices. For (f(U)) we find 

(f)coE = ~U N + ir 1 (trAtiB + tr A T B), 

(/)cse = (2N - l)" 1 (2tr Ati B - tr A R B). (5.7a) 

Similarly, we find for {g(U)) the final result 

<<7>cse = [(2iV-l)(2iV-3)]- 1 

x (AtrAtrBDtrC -2trAD R B R trC - 2 tr A tr BC R D +tvAD R CB R 
+ tr ADCB + tr AC R D R B + tr ADB R C R - 2 tr AC R tr BD R ) 
- [(2A^(2A^- l)(2A^-3)]" 1 

x (2trAtxBDC - tr AC R D R B R +2tiAtiBD R C -trAC R DB R 

+ 2 tr ADB tr C + 2 tr AD R 5 tr C - tr ADBC R - tr AD R £C R 

+ 2 tr A tr BC R D R - tr ADCi? R - 4 tr A tr BD 11 tr C + 2 tr ADB R tr C 

-trAD R S R C R +2trAC R tr5 J D - tr AD R CB -txAC R DB). (5.7b) 



VI. APPLICATION TO A CHAOTIC CAVITY 

We consider the system shown in Fig. |6], consisting of a chaotic cavity attached to two 
leads, containing tunnel barriers. The MxM scattering matrix S is decomposed into iVj x Nj 
submatrices Sy, 

5 = f 511 512 ) , (6.1) 
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FIG. 6. Chaotic cavity (grey) connected to two leads containing tunnel barriers (black). 



which describe scattering from lead j into lead i (M 
given by the Landauer formula, 



Ni + Nj). The conductance G is 



tr s 19 s 



12^12 



= ti C 1 SC 2 S^, G = 2e 2 /h. 
1 — C\ are defined by (Ci)ij 



1 if i 



(6.2) 

j < Nx and 



Gj Go 

The projection matrices C\ and C 2 
otherwise. 

In the absence of tunnel barriers in the leads, S is distributed according to the circular 
ensemble. The symmetry index (3 6 {1,2,4} distinguishes the COE {(3 = 1), CUE (/3 = 2), 
and CSE {(3 — 4). Calculation of the average and variance of G is straightforward |5|||, 

PN X N 2 



(G/G ) 
v&rG/Go 



PM + 2-P' 

2(3N 1 N 2 {f3N 1 + 2- f3){(3N 2 + 2-/3) 



(6.3) 
(6.4) 



((3M + 2- 2(3){(3M + 2- (3) 2 {(3M + 4 - (3) ' 

In the presence of a tunnel barrier in lead i with reflection matrix r^, the distribution of 
S is given by the Poisson kernel [j37H40|| , 

.(J3M+2-/3) - ( n \ 

r 2 ■ 



P(S) oc det(l - S^S) 



S 



(6.5) 



SS = S 



The sub- unitary matrix S is the ensemble average of S: J dSP(S)S = S. The eigenvalues 
are the transmission eigenvalues of the tunnel barriers. The fluctuating part 
S of S can be decomposed as 

5S = T'(l - UR'^UT, (6.6) 

where T, T", and i?' are M x M matrices such that the 2M x 2M matrix 

(6.7) 




is unitary. The usefulness of the decomposition ( |6.6| ) is that U is distributed according the 

In the presence of time-reversal symmetry, we further have 
and U = U T . Physically, U corresponds to the scattering 



circular ensemble |pllp7l^0 



>/T 



S = S T , V = T T , R! = R r 
matrix of the cavity without the tunnel barriers in the leads and £ corresponds to the 
scattering matrix of the tunnel barriers in the absence of the cavity Hi~9| , |40| . 

The decomposition ( |6.6j ) reduces the problem of averaging S with the Poisson kernel to 
integrating U over the unitary group. Because the conductance G is a rational function of U, 
this average can not be done in closed form for all M. For N±, N 2 ^> 1 a perturbative calcu- 
lation is possible. In this section we will compute the mean and variance of the conductance 
in the large- N limit, using the diagrammatic technique of the previous sections. 
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T * R ^ * R * * T 

<■ ■ -C < • ■ ■ ■ ■ ■ -O • ■ ■ -CL- 

■ ■ ■ ■ ■ -O^* ■ ■ -O • ■ ■ >(f* 

r Rf FT T 



a 



T't * Fft * R ^ * * V 

., <ra:m : 

T Ff FT T 

FIG. 7. Top: Diagrammatic representation of the function f n (U) in Eq. ( |6.9[ ); Bottom: Ladder 
diagram with the largest number of cycles, which gives the 0{N) contribution to the average 
conductance. The arrows are omitted if the direction of the diagram is not ambiguous. 

A. Average conductance 

According to the Landauer formula (|6.2| ) the average conductance is given by 

(G/G ) = (tr C 1 8SC 2 5S^, (6.8) 

where we have used that (5S) = 0. Expansion of the denominator in the decomposition 
(|6.6|) of 6S yields the series 



(G/G ) = J2(UU)), (6.9a) 

n=l 

f n (U) = tr C l T\UB!) n - l UTC 2 T ] U ] (R'^U^^T'l (6.9b) 



The average of the polynomial function f n (U) can be calculated diagrammatically. We 
represent f n (U) by the top diagram in Fig. 0. The average over the matrix U is done as 
follows. 

The leading contribution, which is of order M, comes from the diagrams with the largest 
number of T- and [/-cycles. For a polynomial of the type flS.8| ) (all U's are on one side of the 
W's), these diagrams have a "ladder" structure (see bottom diagram in Fig. ^). The ladder 
diagrams contain n [/-cycles and n + 1 T-cycles. Their weight is WJ 1 = M~ n + (9(M _n_1 ), 
resulting in 

(/„([/)) = M- n trT ,t CiT'(tri? , i? ,t )"- 1 trTC 2 T t + C(l). (6.10) 

Summation of the series fl6.9| ) yields (G) to leading order in M, 

_ (trrt gi r)(trTC 2 Tt) 

{G/Go) ~ M-trR'W + 0{l) 

= (^-trffl^-tr^ +0(1) (6n) 

In the second equality we have used the unitarity of the matrix E defined in Eq. fl6.7p. 
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R't 



R't 



•7^ : " " ™ : " 



FIG. 8. Top and middle: Two maximally crossed diagrams contributing to the 

weak-localization correction to the average conductance. The right and left parts of the diagram 
have a ladder structure; Bottom: The maximally crossed part of the top diagram redrawn as a 
ladder diagram. 



The weak-localization correction is the 0(1) contribution to (G). In general, an 0(1) 
contribution to the average conductance can have two sources: (i) a higher order contribution 
to the weight W cl ,..., Ck of the leading-order diagrams, and (ii) higher order diagrams. In the 
absence of time-reversal symmetry both contributions are absent: (i) W\ = M _1 has no 
0(M~ 2 ) term, and (ii) there are no diagrams of order 1. 

The situation is different in the presence of time-reversal symmetry. We discuss the case 
/3 — 1 in which there is no spin-orbit scattering. The case /3 — 4 then follows from the 
translation rule of Sec. [V|. In the presence of time-reversal symmetry, (i) the coefficient 
Wi = M~ l - M~ 2 + ... has an 0(M~ 2 ) term, and (ii) there are diagrams of order 1. The 
first contribution is a correction nM~ n ~ l to the weight M~ n in Eq. ( |6.10| ). Summation over 
n yields the first correction to Eq. (|6.11|) , 



(trT't^rXtrT^Tt) 
6Gl = (M-tri^'t)2 • (6 - 12) 

The second contribution is from diagrams which are obtained from the ladder diagrams by 
reversing the order of the contractions in a part of the diagram. The central part of the 
diagram is "maximally crossed" , the left and right ends are ladders (see Fig. ||). In disordered 
systems, the ladder diagrams are known as diffusons, while the maximally crossed diagrams 
are known as cooperons. The maximally crossed diagrams are not allowed in the absence of 
time-reversal symmetry, because dots of different color are connected by thin lines (violating 



rule 2 in Sec. [TTI|) . A maximally crossed diagram can be redrawn as a ladder diagram by 
nipping one of the horizontal lines along a vertical axis (bottom diagram in Fig. |8|). 

In the maximally crossed diagrams all cycles but one have minimum length. The cycle 
with the exceptional length can be a [/-cycle (top diagram in Fig. H), or a T-cycle (middle 
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diagram). To evaluate these diagrams, we need to introduce some more notation (see Fig. 
H). We denote the left and right ladder diagrams by matrices F L and F R , 

oo 

F L = T'^dT + Y M- n (trT ,t C 1 T')(tr J R ,t J R , ) n ~ 1 ^ /t ^ / 

71=1 
OO 

F R = TC 2 T ] + Y M- n E!R!\tiR!B! ] ) n -\tiTC 2 T ] ) 

n=l 

The scalars fuu an d Jtt represent the maximally crossed part of the diagram, 

f TT = £ M- n (tiR'R^) n+1 = a f trjR 5L (6.14a) 

OO 1 

fuu = Y M~ n ~ 1 (\xR'R^) n = — — -. (6.14b) 

We used the symmetry of R' to replace R' T by R'. With this notation we may draw the 
contribution 5G 2 to the weak-localization correction from the maximally crossed diagrams 
as in Fig. |10[ It evaluates to 

SG 2 = -M- 3 tiFJ TT tiF R + tTFJ uu Fl (6.15) 

The total weak-localization correction 5G = 8G\ + 5G 2 becomes 

SG = -(tr T f T)~ 3 [(tr C 2 T f T) 2 trCi^T) 2 + (tr dT ] T) 2 tr C 2 (T f T) 2 ] . (6.16) 

Since T^T = 1 — SW has eigenvalues r„, we may write the final result for the average 
conductance in the form 



<G/G„>- 9191 



9i + 9'i \ P, 

Nj M 

9 P = YK, 9' P = Y K- (6-18) 

n=l n=l+JVi 

(The (3 = 4 result follows from the translation rule of Sec. |V|.) The first term in Eq. 
06.17] ) is the series conductance of the two tunnel conductances G gi an d G g[. The term 
proportional to 1 — 2/(3 is the weak- localization correction. In the absence of tunnel barriers 
one has g p = iV\, g' = N 2 , and the large-M limit of Eq. ( |6\3| ) is recovered. In the case of two 
identical tunnel barriers (Nx = N 2 = M/2 = N,T n = T n+N for j = 1, . . . , N), Eq. (|6U7D 
simplifies to 



(G/Go) 



l 9i+ i!-^t +0{M ' l) - (619) 

Eq. ( |6.19| ) was previously obtained by Iida, Weidenmiiller and Zuk ||15|| . If all r n 's are equal 
to T, Eq. ( ISlSl) simplifies further to (G/G ) = ±NT + |(1 - 2//3)T. 
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f_ = 



Q- 



c- 



c cc:c: 



<— f 9** o4« ■*■(>«• c>4« ■*■(>#• ■*■(>#• 



UU" 



* 

I 



* * . * * 



9*4t" "9^f ■ "9 



FIG. 9. Diagrammatic representation of Eqs. ( |6.13 ) and ( 6.14| ) 




FIG. 10. Diagrammatic representation of the weak-localization correction 5G2 from the max- 
imally crossed diagrams. The total correction 5G = 5G\ + 5G2 contains also a contribution 5G\ 
from the weight factors [Eq. ( |6.12| )] 
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f kl f mn = 



■ ■c ) •■■■ c) l "-0 •■ 

T 4 . * R'J * R'j * * Tj 

1 ^« ■ ■ ■ ■ ■ -o^* ■ ■ -o • ■ ■ ■ V' 

T R' R' T 

■ ■ •o4+ ■ ■ ■ ■ ■ -o • ■ ■ ■( 

T t R -T R -t T t 

FIG. 11. Diagrammatic representation of a term contributing to G 2 , and hence to the variance 
I2C|) of the conductance. 



B. Conductance fluctuations 



We seek the effect of tunnel barriers on the variance of the conductance, varG = (G 2 ) — 
(G) 2 . We consider (3 = 1 and 2 first, and translate to (3 = 4 in the end. Using the 
decomposition ( |6.6|) we write the variance in the form 



var Gj G = var (tr C 1 5SC 2 SS j ) = covar (foh fn 

k,l,m,n~>l 

fki = tTC 1 T'(UR') k ~ 1 UTC 2 T^U\R^U^ l - 1 T^. 



(6.20a) 
(6.20b) 



Since the number C/'s and ?7*'s must be equal for a non-zero average, covar (f kl , f mn ) = 
(fkifmn) ~ (fki)(fmn) = unless k + in = I + 7i. Diagrammatically, we represent f k if mn by 
Fig. |TT[ The diagram consists of an inner loop, corresponding to fki, and an outer loop, 
corresponding to f mn . The covariance of fki and f mn is given by the connected diagrams. 
We call a diagram "connected" if (i) the partition of the ?7-cycles contains a group which 
consists of ?7-cycles from the inner and the outer part, or (ii) the diagram contains a cycle 
(a [/-cycle or a T-cycle) connecting the inner and outer loops. 

We first compute the contribution from diagrams which are connected only because of 
(i), i.e. diagrams in which all [/-cycles and T-cycles belong either to the inner or outer loop. 
The contribution from such a diagram is maximal, if the [/-cycles are partitioned in groups 
which are as small as possible. The optimal partition consists of groups of size 1, except for 
a single group of size 2, which contains one [/-cycle from the inner and one from the outer 
loop. Furthermore, the total number of cycles is maximal if both the inner and outer loops 
are ladder diagrams. This requires k = I and m = n. The covariance from this diagram is 



covariance = kmSkiS rnn W 1 1 Wj m 2 

x (tr T'^dT) 2 (tri?'i?t) fc+m - 2 (tr TC 2 T^) 2 + 0(M^ V ) 

Summing over k and m we obtain the first contribution to varG/G , 

variance = M -4 (tr i^tr Fr) 2 . 



(6.21) 



(6.22) 



The second contribution, consisting of diagrams in which the inner and outer loops are 
connected by T- or [/-cycles, is of maximal order if the partition of the [/-cycles involves 
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c— 






c— 






c— 









(o) (p) 

FIG. 12. The 16 connected diagrams which contribute to the variance of the conductance. The 
shaded parts are defined in Figs. |9] and [l3|. These diagrams contribute for (3 = 1 and 2. For (3=1 
there are 16 more diagrams, obtained by flipping the inner loop around a vertical axis (diagram 
a-h) or around a horizontal axis (i-p), so that ladders become maximally crossed. 
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diagram 


0=1,2 


= 1 


a 
b 
c 
d 


Wi(trF L ) 2 ft T (trF^ 
W-s (tr F L ) Z }tt (tr F R y 
W 2 (tiF h ) 2 f 2 u tvFl 
tr Hrffuu 


W 2 2 tr Fr tr .Rl /| t tr F L tr Fr 

TT7" 1 D 1 D f J.77IJ.771 

W3 tr F R trF L f TT tr F L tr Fr 
W 2 trF R trF L f 2 u tiF^F R 
tr H T H^ fjjjj 


e 

r 
I 

g 
h 


W 2 trFlfl T (trF R ) 2 
w 3 (tr .tl) Jtt (tr F R ) Z 
trF^trF 2 
tr Hi Hfuu 


W 2 tr F R F\,fu T tr Fl tr Fr 

TJ7" J. T7 i 7? X j 771 j 771 

M/ 3 tr Fr trF L / TT tr F L tr Fr 
t r -^r-^l fuu tr FjJFr 
tr H*H fjjjj 


i 

i 
J 

k 
1 


Wi tr F L tr F R / 2 t tr F L tr F R 
WZ tr Ft tr Fu f%™ tr Ft tr Fr 
W 2 tr Ffl't /rtz/t/r tr F L tr F R 
W 2 tr F L tr F R /ri//i/r tr R' 


W% tr F R tr F L Stt tr F L tr Fr 
Wo tr Ft tr Fr fi 2 ™ tr Fr tr Ft 
W 2 tr TYTfl't /ttj/c/t tr F L tr Fr 
W 2 tr F L tr F R f TU fur tr 


m 
n 
o 
P 


W 2 tr F L tr Fr frufuT tr #fft 
W 2 tr F'ttf / Tl //t/T tr F L tr Fr 
tri?^4 trF/tft 


W 2 tr F R tr F L frufuT tr R fr W 
W 2 tr /riz/t/T tr F L tr Fr 
txH^R'l f^trR^H^ 
tr R^Hffru trH*R' 




TABLE I. Contribution to varG/Go from 


the connected diagrams of Fig. 12. 



TU 







O O 



* * 



n"? '^■"t^"? T^f" ,( ?"^ , "r^" 



o 



o 



* * 



* * * 



f UT = 



□ ■ i-in-i.n:n-in:n:n 



H 



H 






FIG. 13. Diagrammatic representation of Eq. (6.23). 
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only groups of size 1. For (3 = 2 there are 16 connected diagrams of maximal order. They 
are shown in Fig. 12, and their contribution to vsjG/Gq is tabulated in Table |. The shaded 



areas indicate ladder parts of the diagram (see Figs. |9] and |T^). The matrices F L and F R , 
and the scalars fuu an d Jtt are defined in Eqs. (|6.13|) and (|6.14|) . The definitions of the 
matrix H and of the scalars fur and J'tu are 

OO i D/ D/f 

f UT = f TU = V M- n (trR'R'^ n = -r „ — , (6.23a) 
jui jiu v ) M-\xR'RV v ; 

HM~ l (tr F R ) R'T'^CiT' + M" 1 (tr F h )TC 2 T^R' + M~ 2 (tr F h ) {tr F R )R'R'^ R' . (6.23b) 

In the presence of time-reversal symmetry (/3 = 1), the matrix U is symmetric. Diagrammat- 
ically, this means that no distinction is made between black and white dots. In addition to 
the 16 diffuson-like diagrams of Fig. IT2|, 16 more cooperon-like diagrams contribute. These 



are obtained from the diagrams of Fig. 12 by flipping the inner loop around a vertical (Fig. 
|r2]a-h) or horizontal (Fig. IT2i-p) axis, so that segments with a ladder structure become 
maximally crossed. Their contributions are listed in Table |. The contributions from the 
individual diffuson-like and cooperon-like diagrams are different. The total contribution to 
var G from diffuson-like and cooperon-like diagrams is the same. 
The final result for the variance of G is 

var G/G = 2/T 1 (g x + g[) _6 (2g$g? + Agfgf - Ag^gf + 2g 2 l9 [ 4 - 2g x g 2 g't 
+ Zglg't - 2g 1 g 3 g' 1 i + 2g 2 g[ 5 - 2g 3 g[ 5 + 2g\g' 2 - 2g\g' x g' 2 

- ±glg? g'2 + Gg 2 i929?92 + 3^ 2 - ^h's - iglg'A) ■ (6-24) 



One verifies that the large- N limit of Eq. ( |6.4| ) is recovered in the absence of tunnel barriers. 
For the special case of identical tunnel barriers (g p = g' ) , this simplifies to 

var G/G = (8/3 gl)' 1 (2g\ - 2 9l g 2 + 3gj - 2 9l g 3 ) , (6.25) 



in agreement with Ref. ||15|| . If all transmission eigenvalues r n = T are equal, one has 
varG/Go = (8/3) _1 [l + (1 — T) 2 ]. A high tunnel barrier (r <^ 1) thus doubles the variance. 



C. Density of transmission eigenvalues 

The transmission eigenvalues T n e [0, 1] are the N x eigenvalues of the matrix product 
s u s\ 2 . Without loss of generality we may assume that N x < N 2 . The matrix product s 2 iS2i 
then has the same N x eigenvalues as s 12 s 12 , plus N 2 — N x eigenvalues equal to zero. The A^i 
non-zero transmission eigenvalues appear as the diagonal elements of the diagonal matrix T 
in the polar decomposition of the scattering matrix 



( Sll 


S\2 ^ 


-{I 


° 


1 S 2 1 


S22 ) 




W 1 




S= 11 = I t . (6.26) 



Here v and v' (w and w') are N\ x N x (N 2 X N 2 ) unitary matrices and 11 is the N 2 — N\ 
dimensional unit matrix. If N x = N 2 , Eq. ( |6.26| ) simplifies to Eq. ( |1.2| ). 
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F(z) C 1 C 1 S C 2 S T C 1 



1 -2 * 

z — + z -C— ■<)■ ■•— + 



-3 ^ c— * > > C"K3 * j 

c 1 s c 2 s T c 1 s c 2 s T c 1 



z" 1 — z" 2 — 3 * > > ■ .O- + 



F'(z) C 2 C 2 S f C 1 S C 2 

£3 * ^ ^ C— * * > q + 

c 2 s f c 1 s c 2 s f c 1 s c 2 

FIG. 14. Diagrammatic representation of the Green functions for the density of transmission 
eigenvalues. 

Sofar we have only studied the conductance G = Go ^2 n T n . The leading contribution to 
the average conductance comes from ladder diagrams. If we wish to average transport prop- 
erties of the form A = J2 n a {T n ) (so-called linear statistics on the transmission eigenvalues), 
we need to know the density p(T) of the transmission eigenvalues T n . The leading-order 
contribution to the transmission-eigenvalue density is given by a larger class of diagrams, as 
we now discuss. 

The density p{T) = (J2n=ifi(T ~ ^n)) °f the transmission eigenvalues follows from the 
matrix Green function F(z): 

F{z) = (C x {z - SC 2 S^C 1 )~ l ) , (6.27a) 
p(T) = -7r _1 ImtrF(r + ie), (6.27b) 

where e is a positive infinitesimal. We first compute p(T) in the absence of tunnel barriers, 
when the result is known from other methods |4]-|6 ,f4l"|l . Then we include the tunnel barriers, 
when the result is not known. 

In the absence of tunnel barriers, the scattering matrix S is distributed according to 
the circular ensemble, so that averaging amounts to integrating over the unitary group. We 
compute F(z) as an expansion in powers of 1/z, 

oo 

F(z) = J2(Ciz- l (SC 2 S^C 1 z- 1 ) n ). (6.28) 

n=0 

We will also need the Green function 

oo 

F'{z) = (C 2 {z - S^C^)- 1 ) = Y J {C 2 z-\S^C l SC 2 z- 1 ) n ). (6.29) 

n=0 

The two Green functions F and F' are represented diagrammatically in Fig. [14|. A diagram 
contributes to leading order [which is 0(1)] if the number of T-and [/-cycles is maximal. 
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F(z) C 1 C 1 L(z) * F(z) 



E(z) F'(z) 



Z 2 W ic^.U^.Uz^ + ... 

zF'(z) * F(z) zF'(z) 

FIG. 15. Diagrammatic representation of the Dyson equation ( |6.30| ) for F(z). 

That is the case if the diagram is planar, meaning that the thin lines do not cross. The ladder 
diagrams are a subset of the planar diagrams. Planar diagrams have been studied in the 
context of the diagrammatic evaluation of integrals over Hermitian matrices, in particular 
for the Gaussian ensemble fl2|]l7|| . For the Gaussian ensemble, only planar diagrams with 
[/-cycles of unit length have to be taken into account. Summation over all these diagrams 
results in a self-consistency or Dyson equation for F(z), which solves the problem [17]. For 
an integral of unitary matrices, [/-cycles of arbitrary length need to be taken into account, 
as is shown diagrammatically in Fig. |T3[ The corresponding Dyson equation is 

oo 

F(z) = z^d + z- l CJl(z)F(z), = W n [ztiF'(z)] n [trF(z)] n -\ (6.30a) 

n=l 

oo 

F\z) = z~ x C 2 + z~ x C 2 Yl{z)F\z), E'(z) = £ W n [z tr F(z)] n [tr F'{z)] n -\ (6.30b) 

n=l 

In terms of the generating function 

OO 1 

h(z) = W n z n ' 1 = — (VM 2 + Az - hi) , (6.31) 



2^ 



we may rewrite Eq. ( |6.30| ) as 

F{z) = C 1 {z-Y,{z)C i y 1 , E(z) = h(ztrF(z)trF'(z))ztrF r (z), (6.32a) 
F\z) = C 2 (z-Z(z)C 2 y\ E'(z) = h{ztrF{z)%rF'{z))ztxF{z). (6.32b) 



In the derivation of Eq. ( |6.32| ) we did not use the particular form of the matrices C\ and 
C 2 . As a check we may choose C\ = C 2 = 1, so that F(z) = F'(z) = [z — and verify 
that Eq. QOg ) holds. 

The solution of Eq. (|6.32|) is 

tpn Ni ~ N 2 \jM 2 z - ( Nl - N 2 y 

tr F(z) = + ^_ , (6.33a) 



^ Z) ^ ^ M \-^ . (,33b) 

IZ 2Zy/z — 1 
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The resulting density of transmission eigenvalues is 



p(T) = lnTVl _7 g(r ~ rmm) ' Tmin= "'V ' (04) 

in agreement with Refs. |B|)[6|j41| . (The function 6(x) — 1 if x > and if x < 0.) 

The weak-localization correction to p(T) follows from the 0(M~ 1 ) term in the large-M 



\2 



expansion of F(z). As in Sec. [VI A| , it has two contributions: which is due to the sub- 

leading order term in the large-M expansion of W n , and 5F 2 (z), which is due to diagrams of 
order 0(M~ l ). In the absence of time-reversal symmetry, both contributions are absent. In 
the presence of time-reversal symmetry, the sub-leading order term SW n = — M~ 2n (— 4) n ~ 1 
in the large-M expansion of W n [cf. Eq. ( f4.4| )] yields a sub-leading order contribution 5h to 
the generating function h, 

oo 

Sh(z) = J2 SWnZ* 1 - 1 = -(M 2 + 4z)-\ (6.35) 

n=l 

from which we obtain 

tv5F 1 (z) = \{z - T min ) _1 - \{z - I)' 1 . (6.36) 

The contribution 5F 2 (z) comes from diagrams in which thin lines connect black and white 
dots. Each such diagram contains the product C\C^ which vanishes. Hence, the 0(M~ l ) 
contribution to F(z) consists of 8F\{z) only. The resulting weak- localization correction to 
the transmission eigenvalue density is 

6p ^ = ^/ ~ Tmin -e)-S(T-l + e)] , (6.37) 

in agreement with Refs. Jl|||. 

We now include tunnel barriers in the leads. Motivated by Nazarov's calculation of the 
density of transmission eigenvalues in a disordered metal ^2j, we introduce the 2M x 2M 
matrices 

C.[^).m-(^). (,3Sa) 
T'=^' T °„), R'=^ (6.38b) 
Analogous to Eq. (|6.6|) , we decompose S = S + SS, where S = (S) and 

5S = T'(l - UR') _1 UT, U = ( U A J (6.39) 




is given in terms of a matrix U which is distributed according to the circular ensemble. 
Because S, Ci, and C 2 commute and C\C 2 = 0, we may replace S by SS in the expression 
( |6.27a| ) for F(z). The result for the matrix Green function F(z) is 
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F(z) = (2Z)- 1 £ (C ± CT' [l - U(R' ± TCTV 1/2 )1 UTC^ 1/2 
= (2Z)- 1 J] [C ± A±(F± - X±)B±] . 



(6.40) 



In the second equation we abbreviated X± = R' ± TCT'r 1 / 2 , F± = (X±(l - UXi)" 1 ), 
and defined A± and B ± such that A ± X ± = CT', X ± B ± = TCz~ 1/2 . 

After these algebraic manipulations we are ready to compute F± by expanding in planar 
diagrams. The result is a Dyson equation similar to Eq. ( |6.30|) , 



F± = X± (1 + S±F±) , S± = £ W n (VF ± 



,2n-l 



(6.41) 



n=l 



where the projection operator V acts on a 2M x 2M matrix A as 



A 



A-21 A 2 2 



VA 



TL AI tiA 12 
U M tr^ 21 



(6.42) 



11m being the M x M unit matrix. The presence of the projection operator V in Eq. ( |6.41| ) 
ensures that the planar diagrams contain only contractions between U (the 1, 1 block of U) 
and (the 2, 2 block of U). In terms of the generating function h we obtain the result 



F = (2z)- 1 (C ± CT'(1 - S ± X ± )- 1 S ± TC^ 1/2 ) , 
S± = (PX±(1 - SA)- 1 ) /i ((7>X±(1 - SiXi)- 1 



(6.43) 
(6.44) 

It remains to solve the 2x2 matrix equation (|6.44[) . We could not do this analytically 
for arbitrary Tj, but only for the case of two identical tunnel barriers: Ni = N 2 = \M = N, 
Tj = Tj +N (J = 1, 2, . . . , N). The solution of Eq. (|6.44| ) in that case is 



±(Vz- 



1m 
1m 



independent of the T^'s. The trace of the Green function is 

( ) = f 2(1 - r,)(^ - y^T) + Tj/y/z^l 

y) fr{ 2z(i-r i )( v ^- v ^T) + r J Vi 

and the corresponding density of transmission eigenvalues is 



(6.45) 



(6.46) 



r,-(2 - r,; 



N 

p(T) = V 

J= i 7r(r 2 - 4TjT + 4T)JT(1 - T) 



(6.47) 



As a check, we note that p(T) -> N5{T) if Tj -» for all j, and p(T) -> iW _1 [T(l — T)] 
if Tj — > 1 for all j [in agreement with Eq. fl6.34 )]. 



-1/2 
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FIG. 16. Conductor consisting of a normal metal (grey) coupled to one normal-metal reservoir 
(N) and one superconducting reservoir (S). The conductor may consist of a disordered segment or 
of a quantum dot. 

VII. APPLICATION TO A NORMAL-METAL SUPERCONDUCTOR JUNCTION 

As an altogether different application of the diagrammatic technique, we consider a junc- 
tion between a normal metal (N) and a superconductor (S) (see Fig. [L6|). At temperatures 
and voltages below the excitation gap A in S, conduction takes place via the mechanism of 



Andreev reflection [43]: An electron coming from N with an energy e (relative to the Fermi 
energy Ep) is reflected at the NS interface as a hole with energy — e. The missing charge of 
2e is absorbed by the superconducting condensate. We calculate the average and variance 
of the conductance, for the two cases that the NS junction consists of a disordered wire or 
of a chaotic cavity. 

Starting point of the calculation is the relationship between the differential conductance 
G^s{eV) = dl/dV of the NS junction and the transmission and reflection matrices of the 
normal region ||44|| , 

G NS (e) = ^tr (?(£) [1 + r'(-eyr'(e)}- 1 t(-ey) (t'(e) [1 + r'(-e) V^f 1 f(- £ )*) f . (7.1) 

This formula requires eV C A <C Ep and zero temperature. The reflection and transmission 
matrices are N x N matrices, which together constitute the 2N x 2N scattering matrix S. 
Using the polar decomposition ( |1.2| ) we may rewrite the conductance formula ( |7.1| ) as 



,2 



4e^ 

G m (e) = — tr 
h 



T + (l + U + y/l ~ T_U*_y/l ~ T+) 1 U + 

(7.2) 



• T_u\ ( 1 + J I -T + u T Jl -T_u\ ^ 



where T± = T(±e) and u± = w'(±e)w(^fe)* . In the presence of spin-orbit scattering, S is a 
matrix of quaternions, and the transpose should be replaced by the dual. In what follows, 
we will consider the case of no spin-orbit scattering. Spin-orbit scattering (considered by 
Slevin, Pichard, and Mello |[45|| ) will be included at the end by means of the translation rule 
of Sec. 0. 

Averages are computed in two steps: first over the unitary matrix u, then over the matrix 
of transmission eigenvalues T. Four cases can be distinguished, depending on the magnitude 
of the magnetic field B and voltage V relative to the characteristic field B c for breaking time- 
reversal symmetry (T) and characteristic voltage E c / e for breaking electron-hole degeneracy 
(V) 
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u+ VR u. VR + u + m+ 



■o 
■o 



u+ VR. u. VR + u + 

FIG. 17. Ladder diagram for the 0{N) contribution to (Gns)- We defined R± = 1 — T±. 

0. eV <C E c , B <C B c -t==> T and V are both present: Then u± may be approximated by 
the unit matrix, so that one only needs to average over the transmission eigenvalues. 
This case has been studied extensively |47] and does not concern us here. 



1. eV <C E c , B ^> B c T> is present, but T is broken: Then we may neglect the 
e-dependence of S, so that u + = w_ = u. According to the isotropy assumption, u is 
uniformly distributed in U(N). 

2. eV ^> E c , B ^ B c T is present, but T> is broken: Then we may consider S(e) and 



S(—e) as independent unitary symmetric matrices. Hence u_ 
distributed in U(N). 



ill 



u is uniformly 



3. eV ^> E c , B ^> B c -<=>- both T and T> are broken: Then u + and U- are independent, 
both uniformly distributed in U(N). 

We compute the average and variance of the conductance for cases 1,2, and 3. 



A. Average conductance 

We start with the computation of the average conductance (Gns)- We first perform the 
average (• ■ •}„ over u± and then over T±. To leading order only ladder diagrams contribute, 



see Fig. 17. The result is the same for cases 1,2, and 3: 



(G* N s/Go) u = 2N 



Wi- 



th- + 7i_ - T-\+T^_ 



+ 0(1), 



l+'l- 



N 



^trT* = -£?-(+,). 



N 



(7.3a) 
(7.3b) 



The 0{1) contribution 5Gns is different for the three cases. 

Case 1, absence of T and presence of V. We put u± = u, Tk± = t^. For normal metals, 
the 0(1) contribution 5G to (G) vanishes if T is broken. However, in the NS junction an 
(9(1) contribution remains [p6fl . The diagrams which contribute to 5Gns have a maximally 
crossed central part, with contractions between [7's and C/*'s on the same side of the diagram 
(Fig. [Tj|, top). The left and right ends have a ladder structure. In the Hamiltonian approach, 
a similar maximally crossed diagram has been studied by Altland and Zirnbauer ^7|, who 



call it a "symplecton" . In total four diagrams contribute to <5Gns> see Fig. [L9[ The building 
blocks of the diagram have the algebraic expressions 
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Vr, 



Vr 




Vr + u + Vr. u. Vr + u + Vr. 



Vr 



■+ 



Vr. 



«■■■■» * 



Vr, 



Vr. 



Vr. 



FIG. 18. Maximally crossed diagram for the 0(1) correction to (Gns) i n the absence of 
time-reversal symmetry and presence of electron-hole degeneracy (top). The right and left parts 
of the diagram have a ladder structure. The central part may be redrawn as a ladder diagram 
(bottom). 



F ± = T± + (1 - T±) tr T± tr (1 - T T ) ]T N~ 2j ~ 2 [tr (1 - T+) tr (1 - T_)] j 

3=0 

= {n± + T±T lT - T 1+ Ti_) (n + + 7i_ - Ti + Ti_) _1 , 

oo 

-(1 - T T ) trT± £ N~ 2j - 1 [tr (1 - T+) tr (1 - T_)] j 
- (ti± - n ± T T ) (n + + 7i_ - ri + Ti_) _1 , 



/i 



# ± = iN- l T±^l - T± trF T - iiV~ 2 (l - T±)^/l - T± tr F T tr F£, 
/tt± 



-tr (1 - T±) J] ^ ^ (1 - r+) tr (1 - T_)} 3 

j=0 

-N(l - Ti±) (n + + Ti_ - Ti+Ti-) -1 , 

oo 

-tr (1 - T±) J] N~ 2 ^ 2 [tr (1 - T+) tr (1 - T_)] j 

j=0 



N \i- r 1± ) [n+ + n_ - r 1+ n_ 



-i 



,-i 



(7.4a) 



(7.4b) 
(7.4c) 



(7.4d) 

(7.4e) 
(7.4f) 



Capital letters indicate matrices, lower-case letters indicate scalars. The subscripts ± are 
omitted from Fig. [L^ because of electron-hole degeneracy. The 0{1) correction <5Gns repre- 
sented in Fig. O equals 
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TT 

FIG. 19. Diagrams for the 0(1) correction to (Gns) i n the absence of time-reversal symmetry 
and presence of electron-hole degeneracy. 



8G NS /G = 8flj V tTiHVT=T + 4W 2 f TT [{trF) 2 + (trF') 



\t x - 4rf + 4rf 



5Tn 



(7.5) 



We still have to average over the transmission eigenvalues. We use that the sample-to- 
sample fluctuations r k — (r^) are an order 1/N smaller than the average. (This is a general 
property of a linear statistics, i.e. of quantities of the form A = J2n a (T n ), see Ref. Q.) 
Hence 



(f(r k )) = f((r k ))[l + 0(N-% 



(7.6) 



which implies that we may replace the average of the rational functions (|7.3|) and (|7.5|) of 
the Tfc's by the rational functions of the average (r^). This average has the 1/N expansion 



(r k ) = (r k ) + O(N- 2 ), 



(7.7) 



where (rfe)o is O(N ). There is no term of order N 1 in the absence of T. The average over 
T of Eqs. ( |7.3|) and ( |7.5|) becomes 

(G m /G ) = _ 8<n>o- 4< ri >g + 4<^>g -8<7- 3 > + x 



(n)o 



(n>o(2 - <rOo) a 



Case 2, presence of T and absence of P. We put u + 
comes from the maximally crossed diagrams of Fig. |2^, 



u_ 



u. The 0(1) correction 



SG NS /G = 2W 2 tr F + f TT _ tr F'_ + 2 W 2 tr Fj. / Tr+ tr F_ 
+ 2trF + / I/[/ _Fi T + 2tvF[f uu+ F^. 



(7.9) 



Averaging over the transmission eigenvalues amounts to replacing r k ± by its average, r k ± — > 
(Tk)o + N^St}. + 0(N~ 2 ). (The average of r k ± is the same for +e and — e.) Because T is not 
broken there is a term of 0(iV -1 ) in this expression. We find for the average conductance 
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FIG. 20. Diagrams for the 0(1) correction to (Gns) m the absence of electron-hole degeneracy 
and presence of time-reversal symmetry. 



(Gns/Cq) 



2iV(7i)o 

2 - (n)o 



+ 



4(ri>g-4(r 1 )g-4(r 2 > + 4(n) (r 2 > 
<n>o (2 - (rOo) 3 



+ C(A^- 1 ). (7.10) 



Case 3, both T and X> broken. Because w + and w_ are independent, there are no diagrams 
which contribute to order 1. The average conductance is obtained by averaging Eq. ( |7.3|) 
over the transmission eigenvalues, 



(Gns/ 'Gq) 



2N(n) 



+ 0(N- 1 ). 



(7.11) 



2 " (7i) 

From the translation rule of Sec. |V| one deduces that in the presence of spin-orbit scat- 
tering, the leading 0(N) term of the average conductance is unchanged, while the 0(1) 
correction is multiplied by —1/2, in agreement with what was found by Slevin, Pichard and 



Mello 45 



The formulas given above apply to any system for which the isotropy assumption holds. 
We discuss two examples: 

(a) A disordered wire (length L, mean free path £, number of transverse modes N), 
connected to a superconductor. We use the results 



(n)o 

( T 2)o 

Sn 



-L/£ 
\(l + £/L)-\ 



'I 



Ki + L/ey' + Hi + L/iy 



(7.12a) 
(7.12b) 
(7.12c) 



We assume f <L<C N£ and neglect terms of order L/N£ and £/L but retain terms of order 
1 and N£ p /LP (p > 1). Substitution of Eq. ( fTT2| into Eqs. (|7J), (|730Q , and Q73ID yields 



(Gns/Gq) 



f ^(l + L/Q- 1 -1 + 4/tt 2 (X>, T), 
I N(l/2 + Lj£)- X - 1/3 (V, no T), 
* N(l/2 + L/e)' 1 - 2/3 (no V, T), 



(7.13) 



[ N(l/2 + Lj£)- X 



(no £>, no T). 
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The result in the presence of both T and D has been taken from Refs. [f!9|,[5(J • 111 the presence 
of spin-orbit scattering, the 0(N) term is unchanged, while the 0(1) term is multiplied by 
-1/2. 

(b) A chaotic cavity without tunnel barriers in the leads. Lead 1 (with Ni modes) is 
connected to a normal metal, lead 2 (with N 2 modes) to a superconductor. An asymmetry 
between N\ and N 2 appears because the dimension of u± in the polar decomposition (|6.26|) is 
N 2 x N 2 . The N 2 x N 2 matrix T± contains the min (Ni, N 2 ) non-zero transmission eigenvalues 
on the diagonal (remaining diagonal elements being zero). We denote N tot = Ni + N 2 and 
Na = (Nf + 6NiN 2 + N 2 ) l l 2 . The averages (ri) and (r 2 )o and the correction 5t% can be 
computed from the density of transmission eigenvalues [Eqs. ( |6.34j ) and ( |6.37| )1. The results 



arc 



Sri = -^NzN'l (n) = NiNg, (t 2 ) = N 1 {N* a -N 1 N 2 )N£. 



(7.14) 



Substitution into Eqs. (fT5|) , ( PH|) , and (|77TT| ) gives 



(G NS /G ) 



( N tot (l 



N tot /N A ) - SN^NlJNj 



(V, T), 



2iV 1 AT 2 /(Ar ot + N 2 ) - 4JSf 1 N 2 N tat /(N tat + N 2 f (V, no T), 
2N l N 2 /(N tot + N 2 ) - AN 2 NlJ(N tot + N 2 f (no V, T), 
[2N l N 2 /(N tot + N 2 ) (noD.noT). 



(7.15) 



The leading order term in Eq. ( |7.15| ) has also been obtained by Argaman and Zee ||29|| . (The 
case Ni = N 2 was given in Ref. 0). 



B. Conductance fluctuations 

To compute the variance of the conductance, we average in two steps: (•••) = ((•■ ~)u)t, 
where (• ■ ■)„ and (■ • -)t are, respectively, the average over the unitary matrices u± and over 
the matrices of transmission eigenvalues T±. It is convenient to add and subtract ((Gns)u)t) 
so that the variance splits up into two parts, 

var G NS = ((G NS )l) T - (<G NS >X + (< G ns>« - (G NS )l) T , (7.16) 

which we evaluate separately. 

The first two terms of Eqs. ( [7.1 6| ) give the variance of (Gns)u over the distribution of 



transmission eigenvalues. We calculated (Gns)u i n Eq. ( [7.3|) . Since (Gns)u is a function of 
the linear statistic t±± only, we know that its fluctuations are an order 1/N smaller than the 
average. This implies that, to leading order in 1/N, 

in \2\ I m \ \ 2 / 9(G m )u\ / d(G NS )u\ 1 \ 

a,a' = ± \ It \ OT lv' I T 

= 8^(2-^.)^ X {J ^ l l <" 7 ' 

We now turn to the third and fourth term of Eq. 1 6| ) . These terms involve the 
variance {G"^ s ) u — {G NS )^ of Gns over U(N) and subsequently an average over the T n 's. 
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FIG. 21. Diagrammatic representation of G^ s . 




The calculation is similar to that of Sec. |VI B| . We represent G^ s by the diagram in Fig. 
|2"T| . The variance with respect to u± is given by the connected diagrams. We distinguish 
between two types of connected diagrams: (i) diagrams in which the inner and the outer 
loop are connected by a T-cycle or by a [/-cycle, and (ii) diagrams in which the partition of 
the [/-cycles involves a group which consists of a [/-cycle from the inner loop and a [/-cycle 
from the outer loop. The diagrams are similar to those of Fig. and are omitted. The 
final result is 

(G 2 m ) u - (G m )l) T =8Gl(2- (r^r 6 (n} 2 

x (4(702 " *(n)l + 9{n)* - A( n )l + 2(71)6 - 4(71)0(72)0 

+ 2(n)2(r 2 ) - 2(7-1)3(72)0 - 2(ri>*(r 2 ) + 6(r 2 )l - 6(ri) (r 2 )2 

+ 3(ri)^(r 2 )2 - 4(ri)o(r 3 )o + 6(r 1 ) 2 (r 3 ) - 2 (n)* (73)0) 
2 (P, no T), 

x < 2 (T, no P), (7.18) 
1 (no P, no T). 



The sum of Eqs. f [7.17|) and (|7.18| ) equals varGNs, according to Eq. ( |7. 1 6|) . 

In the presence of spin-orbit scattering varGws is f° ur times as small, according to the 
translation rule of Sec. [Vj 

We give explicit results for the disordered wire and the chaotic cavity. 

(a) For the disordered wire one has HID] var r i = TE N ~ 2 > ( T k)o = U e / L ) T & T ( k )/ T ( k+ 
|). Substitution into Eqs. ( [7.1 7| ) and ( [7.18j ) yields the variance 



f 16/15 - 48/tt 4 « 0.574 (P, T), 

„ , 8/15 « 0.533 (P, no T), 

varGNs/Go " 8/15 « 0.533 (T, no P), 

[ 4/15 « 0.267 (no P, no T). 



(7.19) 



The result in the presence of both T and P has been taken from Ref . [0,^1 • ^ both P and 
T are present, breaking T (or P) reduces the variance by less than 10% [ p8| , |5"3|j . 

(b) For the chaotic cavity one has var n = 2Nf/f3N^ ot and (r 3 ) = iV 1 (iV'4 t -2iV t 2 ot iV' 1 iV 2 + 



2iVf iVf )/iV t 5 ot [see Eqs. (0) and ( ggg )]. In combination with Eq. (|7TTJ) this gives 
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128iV 2 iVf (A£ t + 2iV 2 iV 2 )(iV 2 t + 4iV 1 iV 2 )- 4 (V, T), 

r : C , 32iV 2 iV 2 t (iV 2 ot - N 1 N 2 )(N tot + N 2 y 6 (V, no T), 

G NS /G - < _ N 1 N 2 )(N tot + N 2 )~* (T, no 2>), (7 - 20) 

16iV|iV 2 ot (iV 2 ot - iV x iV 2 )(iV tot + iV 2 )- 6 (no P, no T). 

If the coupling between the cavity and the normal metal is weak compared to the coupling to 
the superconductor (iV 2 > Aq), one finds varG NS (£>, T)/var G m (D, no T) = 0(N 1 /N 2 ) 2 . 
In this case breaking T greatly enhances the conductance fluctuations. In the opposite 
case, if the couplings are equal (Ni = N 2 ), one finds vaiG^si?), T)/var Gns^, no T) = 
2187/2084 fa 1.07. In this case breaking T has almost no effect on the conductance fluctu- 
ations. 



VIII. SUMMARY 

We developed a diagrammatic technique for the evaluation of integrals of polynomial 
functions of unitary matrices over the unitary group U(N). In the large- N limit the number 
of relevant diagrams is restricted, which allows for the evaluation of integrals over rational 
functions. We also considered integrals of unitary symmetric matrices, by means of a slight 
modification of the diagrammatic rules. A translation rule was given to relate integrals of 
(self-dual) unitary matrices of quaternions to integrals over (symmetric) unitary matrices of 
complex numbers. 

We discussed two applications: A chaotic cavity (quantum dot) with tunnel barriers in 
the leads and a normal-metal-superconductor (NS) junction. In both cases, the conductance 
is a rational function of a unitary matrix. In the large- iV limit the average conductance is 
given by a series of ladder diagrams. The weak-localization correction consists of maximally- 
crossed diagrams. These two types of diagrams are analogous to the diffuson and cooperon 



diagrams in the diagrammatic perturbation theory for disordered systems |]22| , [23|| . We com- 
puted the density of transmission eigenvalues, where the leading order term is given by 
planar diagrams. Resummation of the diagrams leads to a Dyson equation for the Green 
function, similar to that encountered in the theory of integrals over Hermitian matrices 



12,17 



For the NS junction, the 0(1) correction to the average conductance is non-zero in the 
presence of a magnetic field, because of a different type of maximally crossed diagrams. These 
diagrams are suppressed by a sufficiently large voltage to break electron-hole degeneracy. 
The new type of maximally crossed diagrams explains the coexistence of weak localization 
with a magnetic field |26| and the insensitivity of the conductance fluctuations to a magnetic 



field 

This research was supported by the "Nederlandse organisatie voor Wetenschappelijk On- 
derzoek" (NWO) and by the "Stichting voor Fundamenteel Onderzoek der Materie" (FOM). 

APPENDIX A: WEIGHT FACTORS FOR POLYNOMIAL INTEGRALS 
In Tables |ll] - we list the weight factors V Cu ..., Ck and W / Cl) ..., Cfc for n = c\ + . . . + Ck < 5 



for the CUE and the COE. (Tables of V are also given in Refs. pTlp2j for the CUE and 
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n 


A n (CUE) 


A n (COE) 


1 


N 




N + 1 


2 


N(N — 1) 




N(N + 1)(N + 3) 


3 


iV(iV 2 - 1)(JV 2 - 4) 




(JV - l)JV(jV + l)(iV + 3)(JV + 5) 


4 


iV 2 (iV 2 -l)(iV 2 -4)(iV 2 -9) 


(N 


- 2)(iV - 1)JV(JV + 1)(JV + 2)(JV + 3) 








x (iV + 5)(iV + 7) 


5 


N 2 (N 2 - 1)(N 2 - A)(N 2 - 9)(N 2 - 16) 


(N 


- 3)(N - 2)(N - l)N(N + 1)(N + 2) 








x (7V + 3)(iV + 5)(iV + 7)(iV + 9) 



TABLE II. Denominators A n of the coefficients V C1 ... c , for n = ci + . . . + < 5. 



in Ref. |34| for the COE.) The weight factors are rational functions of the dimension iV 
of the unitary matrix. The denominators A n and B n of, respectively, V Cl ,... )Ck and W cl ,..., Ck 
depend only on n. They are tabulated in Tables |I| and |V[ The numerators A n V clt ___ }Ck and 
B n W cl) ... >Ck are tabulated in Tables |IT| and 0. 
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Cl, ... ,c k 


A n V cl ,..., Ck (CUE) 


A n V Cl _ Ck (COE) 


1 


1 


1 


1,1 


TV 


2 + AT 


2 


-1 


-1 


1,1,1 


-2 + iV 2 


2 + 5 N + N 2 


2,1 


—TV 


-3- N 


3 


2 


2 


1,1,1,1 


6-8A^ 2 + 7V 4 


-32 - 8 N + 28 N 2 + 11 N 3 + N 4 


2,1,1 


4N- N 3 


-4.-18N -9N 2 - N 3 


2,2 


6 + N 2 


24+7 N + N 2 


3,1 


-3 + 2N 2 


10 + 12N + 2N 2 


4 


—5 A" 


—11 — 5 N 


1,1,1,1,1 


78 N - 20 iV 3 + iV 5 


128 - 408 N - 84 N 2 + 59 N 3 + 16 N 4 + N 5 


2,1,1,1 


-24 + 14 7V 2 - iV 4 


92 + 38 N - 43 N 2 - UN 3 - N 4 


2,2,1 


-2N + N 3 


56 + 43 A^ + 12 A^ 2 + A^ 3 


3,1,1 


-18N + 2N 3 


-52 + 40N + 22N 2 + 2N 3 


3,2 


-24- 2iV 2 


-88- 18 N - 2N 2 


4,1 


24 - 5 iV 2 


-7-36N -5N 2 


5 


14 N 


38+ UN 



TABLE III. Numerators A n V cli ... jCk of the coefficients Cjt for n = c\ + . . . + < 5. The 
denominators A n are given Table 0. 



n 


B n (CUE) 


S n (COE) 


1 


N 


AT+1 


2 


N 2 {N 2 - 1) 


AT(A^ + l) 2 (AT + 3) 


3 


A^ 3 (A^ 2 - 1)(AT 2 -4) 


(iV - 1)N(N + 1) 3 (AT + 3)(AT + 5) 


4 


N 4 (N 2 - 1) 2 (N 2 - 4)(N 2 - 9) 


(N - 2)(N - 1)N 2 (N + 1) 4 (N + 2) 






x (N + 3) 2 (N + 5)(N + 7) 


5 


^5(^2 _ 1)2(^2 _ 4)^2 _ g^ N 2 _ jgj 


(N - 3)(N - 2)(N - 1)N 2 (N + 1) 5 (N + 2) 






x (N + 3) 2 (N + 5)(A^ + 7)(AT + 9) 



TABLE IV. Denominators B n of the coefficients W Cl ,..., Ck for n = ci + • • • + Ck < 5. 
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Cl, ... ,c k 


BnWc u ...,c k (CUE) 


BnW Cl ,...,c k (COE) 


1 


1 


1 


1 1 

2 


1 

-TV 


2 

-1 - TV 


1,1,1 

2,1 
3 


8 

-4 7V 
2 TV 2 


32 

-8- 87V 
2 + 4 TV + 2 TV 2 


1,1,1,1 
2,1,1 
2,2 
3,1 
4 


-216 + 144 TV 2 
72 iV - 48 TV 3 
-42 7V 2 + 18 7V 4 
-15 7V 2 + 15 7V 4 
5 TV 3 - 5 TV 5 


-1680 + 6720 TV + 6096 TV 2 + 1152 TV 3 
280 - 840 TV - 2136 TV 2 - 1208 TV 3 - 192 TV 4 

-140 - 116 TV + 384 TV 2 + 592 TV 3 + 268 TV 4 + 36 TV 5 
198 TV + 552 TV 2 + 540 TV 3 + 216 TV 4 + 30 TV 5 

-33 TV - 125 TV 2 - 182 TV 3 - 126 TV 4 - 41 TV 5 - 5 TV 6 


1,1,1,1,1 
2,1,1,1 
2,2,1 
3,1,1 
3,2 

4,1 
5 


-13824 + 4224 TV 2 
3456 7V - 1056 TV 3 
-1248 7V 2 + 288 TV 4 
-480 TV 2 + 240 TV 4 
312 7V 3 - 72 7V 5 
56 TV 3 - 56 TV 5 
-14 7V 4 + 14 7V 6 


-483840 + 297984 TV + 407040 TV 2 + 67584 TV 3 
60480 + 23232 TV - 88128 TV 2 - 59328 TV 3 - 8448 TV 4 
-12096 - 21120 TV + 1152 TV 2 + 18432 TV 3 + 9408 TV 4 + 1152 TV 5 
-3024 + 192 TV + 15072 TV 2 + 18432 TV 3 + 7536 TV 4 + 960 TV 5 
1512 + 4152 TV + 2496 TV 2 - 2448 TV 3 - 3480 TV 4 - 1320 TV 5 - 144 TV 6 
-912 TV - 3376 TV 2 - 4768 TV 3 - 3168 TV 4 - 976 TV 5 - 112 TV 6 
114 TV + 536 TV 2 + 1018 TV 3 + 992 TV 4 + 518 TV 5 + 136 TV 6 + 14 TV 7 



TABLE V. Numerators B n W cl Ck of the coefficients W clj _ >Ck for n = c\ + . . . + < 5. The 



denominators B n are given in Table IV 
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